Double optical spring enhancement for gravitational wave detectors 
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Currently planned second-generation gravitational-wave laser interferometers such as Advanced 
LIGO exploit the extensively investigated signal-recycling (SR) technique. Candidate Advanced 
LIGO configurations are usually designed to have two resonances within the detection band, around 
which the sensitivity is enhanced; a stable optical resonance and an unstable optomechanical reso- 
nance - which is upshifted from the pendulum frequency due to the so-called optical-spring effect. 
Alternative to a feedback control system, we propose an all-optical stabilization scheme, in which 
a second optical spring is employed, and the test mass is trapped by a stable ponderomotive po- 
tential well induced by two carrier light fields whose detunings have opposite signs. The double 
optical spring also brings additional flexibility in re-shaping the noise spectral density and optimiz- 
ing toward specific gravitational-wave sources. The presented scheme can be extended easily to a 
multi-optical-spring system that allows further optimization. 

PACS numbers: 04.80.Nn, 03.65.Ta, 42.50.Dv, 42.50.Lc, 95.55.Ym 



cr. 



> ■ 
^ ; 
o\ . 
o ■ 

cn ; 
o 

00 

O 



X 



I. INTRODUCTION 

The array of the four large scale laser interferometers 
(LIGO [H, VIRGO 0, GEO [| and TAMA H) repre- 
sent the first generation of interferometric gravitational- 
wave (GW) detectors {interferometers for short). Next- 
generation interferometers, such as Advanced LIGO 
plan to use the so-called detuned signal-recycling (SR) 
technique — in which an additional mirror is placed be- 
hind the dark port of a Michelson interferometer in or- 
der to modify the optical resonance structure of the in- 
terferometer: depending on the location and the reflec- 
tivity of this signal-recycling mirror, the eigenfrequency 
and quality factor of the optical resonance can be ad- 
justed JpJaS. As shown theoretically by Buonanno and 
Chen [a [13, mi and experimentally by Somiya et al. 
and Miyakawa et al. , detuned signal- recycling makes 
the power inside the interferometer dependent on the mo- 
tion of the mirrors, creating an optical spring. With 
optical power as high as planned for second-generation 
interferometers, detuned signal-recycling interferometers 
are characterized by two resonances within the detec- 
tion band. One resonance is optical in nature, while the 
other one is due to the optical spring: the eigenfrequency 
of the test masses can be shifted from that of the pen- 
dulum up into the detection band. The general principle 
underlying the optical spring effect in signal-recycling in- 
terferometers is identical to that explained by Braginsky, 
Gorodetsky and Khalili [ij] for a single detuned cavity 
(cf. [Hj), which had been employed in their proposal of 
the optical bar detection scheme [l5| . 

One concern that arises with the use of the optical 
spring is that it always causes instability: depending on 
the sign of the detuning, the optical force either brings 



anti-damping, or creates an anti-spring. Buonanno and 
Chen have shown in Ref. [l^l that one can cope with 
the instability by incorporating a linear feedback control 
which ideally would not modify the noise spectral density 
of the GW detector. In practice, however, the need for 
control inside the detection band, can cause undesirable 
complexity in the control system or additional classical 



Here we propose an alternative way to suppress the 
instability, by injecting a second carrier field from the 
bright port (cf. Fig.[T]). We shall assume these carriers 
have different polarizations (as in Ref. [l^), so that there 
is no direct coupling between the two fields (although 
they both directly couple to the mirrors). Differently 
from Ref. , in which the second carrier does not enter 
the arm cavity at all, in this proposal, the second carrier 
resonates in the arm cavity, but is subject to different 
SR detuning and SR mirror reflectivity. Nevertheless, 
similar to Ref. [l6|, a homodyne detection must be per- 
formed at the dark port around each of the two carriers, 
with the two outputs combined with appropriate fllters. 
The main purpose of this second carrier is to create a 
second optical spring that forms a stable optical spring 
together with the first one — even though each individual 
optical spring, acting alone, would be unstable. Such a 
stable, Double Optical Spring (DOS) is possible at least 
in two ways. The first way {weak stabilization) relies on 
the observation that the ratio between the optical spring 
constant's real part (strength of spring/anti-spring) and 
its imaginary part (strength of damping/anti-damping) 
depends on the detuning frequency of the carrier — the 
detunings of the two carriers can be arranged such that 
the first one has a stronger spring and a weaker anti- 
damping, while the second one has a weaker anti-spring 
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FIG. 1: Schematic plot of a power- and signal-recycled 
Michelson interferometer with arm cavities and double optical 
spring, which is realized by injecting two carriers with orthog- 
onal polarizations. Both carriers are split at the beam split- 
ter (BS) and transmitted through the input mirrors (ITM) 
into the arm-cavities consisting of ITM and the end mirrors 
(ETM) ; they both resonate in the arm cavity. Two additional 
cavities are realized by the signal-recycling mirror (SRM) and 
the power- recycling mirror (PRM). 



and a stronger damping. In this case, stabilization can be 
achieved with a weak second carrier, which does not mod- 
ify the sensitivity of the interferometer by much. The sec- 
ond way {annihilation) requires the two carriers to have 
equal power and exactly opposite detunings, such that 
their optical springs exactly cancel each other, and the 
test-masses' differential motion behaves as a free mass 
again. Interestingly, this configuration gives exactly the 
noise spectrum that had been expected by the GW com- 
munity [r^ at the time an interferometer with detuned 
signal-recyc ling was not thoroughly studied, i.e. before 
Refs. i, M, \lM were published. 

Although straightforward to understand, weak stabi- 
lization and annihilation are by no means the only pos- 
sibilities. In fact, an additional benefit of the DOS tech- 
nique is that it increases the flexibility in shaping the 
noise curves: for any specific source, the noise spectrum 
can be optimized correspondingly over the parameter 
subspace of the two carriers, subject to the constraint 
that the resulting dynamics must be stable. In this pa- 
per, we will carry out an optimization for neutron-star- 
neutron-star binary inspirals — using both the current 
noise budget of Advanced LIGO and a plausible noise 
budget for interferometers right beyond the second gener- 
ation. Note that the parameters for the Advanced LIGO 
configuration as well as for the other single optical spring 
configuration have been obtained by using the same op- 
timization. 



The DOS technique can also be used for the stabiliza- 
tion of the optical-spring ponderomotive squeezer, which 
generates frequency independent squeezed vacuum below 
the optical spring frequency [H, [l^ . The stable optome- 
chanical resonance has already been demonstrated exper- 
imentally [20| . 

This paper is organized as follows: In Sec. [Til we shall 
motivate the stabilization process and studying the clas- 
sical dynamics of the double optical spring systems. In 
Sec.|IlT]we recall the necessary basics in order to calcu- 
late the input-output relation of a DOS interferometer. 
In Sec. IIVI different applications of the DOS are discussed 
and our scheme is extended to a a multi-optical-spring 
system. In Sec.|V]we summarize our main conclusions. 



II. CLASSICAL DYNAMICS 

In this section, we consider the classical dynamics of 
double-optical-spring stabilization. For the mechanical 
degrees of freedom, we only consider the differential mode 
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between the interferometer's test-masses, which is sensed 
at the dark port by both carriers. Before coupling to 
the light, the differential mode has an eigenfrequency of 
« IHz (that of the suspension system), and effective 
mass of m/4, where m is the mass of each individual 
mirror. Since the pendulum frequency is far below the 
detection band, we will simply treat the mirrors as free 
masses. For a carrier with angular frequency loq that is 
resonant in the arms, when the mirrors are held fixed, 
the optical resonant frequency of the differential optical 
mode (to be precise, the one that is closest to this carrier) 
is given by — A — ie, where, in terms of interferometer 
parameters, e and A are given by (cf. [ll|) 
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where 70 = Tc/(4L) is the half line width of the arm cav- 
ity (T the power transmissivity of the ITM, c the speed 
of light, and L the arm length), psR is the amplitude re- 
flectivity of the SRM, and (j) the detuning phase of the 
carrier in the SR cavity (single trip). In reality, when 
two carriers are both resonant in the arm cavity, their 
detuning phases in the signal-recycling cavity must differ 
by 
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where (Ai^)FSR is the free spectral range of the arm cavi- 
ties. This constraint must be taken into account in prac- 
tical designs of DOS interferometers. 
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Now with optomechanical coupUng. let us first consider 
a single detuned carrier. Treating the mirrors as free 
masses to start with, the classical equation of motion of 
the differential mode can be written in frequency domain 
as: 
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where m/4 is the reduced mass, Fext is any external clas- 
sical force. The frequency dependent optical spring con- 
stant is given by [ll| 
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where 9 is given by 
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with P the carrier's circulating power in the arms. Note 
that 9 has units of Hz"^. 

Before treating a full-power interferometer, it is in- 
structive to first draw our attention to a weakly coupled 
regime, assuming that optical frequency scales (effective 
detuning A and effective bandwidth e) are much larger 
than the resulting optomechanical resonant frequency. In 
the weakly coupled regime, one can expand the optical 
spring constant (cf. e.g. Ref. [l^l) as 
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= K-inv . (8) 



where K and F are both real constants. Analogous to a 
mechanical spring, K describes the restoring force while 
r denotes the damping. Real and imaginary parts of the 
optical spring constant are proportional to 9, which is 
in turn proportional to the carrier power P. Inserting 
Eq. ([5]) into Eq. stability requires 



K>Q and E > , 



(9) 



or basically positive spring constant and positive damp- 
ing. 

Erom Eq. ([8]) , it is straightforward to deduce that the 
stability condition @ can never be fulfilled by a single 
optical system, since a positive detuning produces always 
anti-damping (F < 0) while a negative detuning always 
comes along with an anti- restoring force {K < 0). In 
Fig-El for fixed circulating power and effective line width 
(e), as the effective detuning (A) shifts from ~oo to oo, 
we plot the trajectory mapped out by {K, F). We use 
two different powers, with outer trajectory correspond- 
ing to a higher power. Indeed, these oo-shaped trajec- 
tories are confined within quadrants with • F < 0. In 
case of a double optical spring, each individual spring, as 
we change its detuning frequency, has its own oo-shaped 
trajectory. When two optical springs combine, their com- 
plex spring constants add up, which can be depicted by 
a vector addition in Eig.[21 By adjusting the detunings 



/ -o.4iV Statically 

A 113 /"■""^ 
\ 1 32\ \ 

\ \ 

1 88 X \ ^ 


Stable 

/ ^ 




\ 188 
\ 1 60 \ 


Anti-Stable 


0.95 1 \ 




0-24V"^ 1.13* 




DynamicallyX^jj 1 




unstable >>^^ 



K 

FIG. 2: Real and imaginary parts of the spring constant Kos 
in the weakly coupled regime [Eq. For each trajectory 

the circulating power and bandwidth are fixed to a certain 
value while detuning varies from — oo to oo. Outer trajectory 
corresponds to a higher circulating power. Example values of 
A/e are marked on the trajectories. 



for first and second carrier it is possible to find many 
stable compositions, with one of them depicted in Fig. [21 
In this configuration, a relatively strong optical spring 
is stabilized by a relatively weak anti-spring generated 
by a lower power. This is possible because the stronger 
optical spring is generated by a carrier with relatively 
high optical quality factor, |A/e|, which tends to yield 
a stronger restoring (or anti-restoring) than damping (or 
anti-damping), while the weak anti-spring is generated by 
a carrier with low optical quality factor, which tends to 
yield a stronger damping (or anti-damping) than restor- 
ing (or anti-restoring). A lower optical power of the sec- 
ond carrier allows the damping of the second spring to 
match that of the first one, while makes the anti-restoring 
of the second spring much weaker than the restoring of 
the first spring. Mathematically, this weak stabilization 
can be summarized as 
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|f(i)| ^ |f(2)| =^ 17^(1)1 » |i^(2)| 



(10) 

As another (rather extreme) example of DOS stabiliza- 
tion, we note that when the two carriers have the same 
power and bandwidth, but the opposite detuning, their 
optical spring constants exactly cancel with each other. 
In fact, this cancelation, or annihilation, is valid for an 
arbitrarily strong coupling, cf. Eq. ([6]). Mathematically, 
this can be summarized as 



Kosi9, e, A; n) + K,si9, e, - A; Q) - 0. 



(11) 



As stated before, the stability condition given in 
Eq. ^ is an approximation only valid in the weakly cou- 
pled regime. A more precise statement regarding the sta- 
bility of the two carrier system is given by the condition 
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that all roots of the characteristic equation 

-!^n' + K^l\n) + KSHn) = Q (12) 

must have negative imaginary parts. In Fig.[3l we explore 
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FIG. 3: Example of DOS stabilization process, as illustrated 
by trajectories of the optomechanical resonant frequency in 
the complex plane [Eq. (|12|) ] . The test masses start off as free 
masses (gray dot). As the first carrier light (e*^^ — 2tt 120Hz 
and A'^' = 2tt 290Hz) is turned on, it causes an upshift of 
the mechanical resonant frequency , as well as a mild anti- 
damping. The trajectory ends at the blue dot, with a circu- 
lating power of 800 W. Subsequently, the second carrier is also 
turned on, bringing a damping while slightly downshifting 
the optomechanical resonant frequency. For e*'^^ = 2-k 5Hz, 
A*-^' = —2tt 55IIz, the trajectory ends at the red dot, with 
a circulating power of 8 kW; while for e'^-* = 2n GOHz, 
A''^' = —2tt 60Hz the trajectory ends at the green dot, with 
circulating power of 80 kW. 

high-power DOS stabilization by tracing the real and 
imaginary parts of the optomechanical eigenfrequency, 
obtained numerically solving Eq. (fT2|) . We first consider 
single-optical-spring configurations with power increas- 
ing from to 800 KW (from left most dot to top dot). 
In this case both the real part of the optomechanical res- 
onant frequency and the anti-damping increases. Then 
we fix the first carrier at 800 kW and increase the sec- 
ond carrier (up to 8kW and 80 kW, respectively, for two 
different choices of linewidth and detuning), which sta- 
bilizes the system by adding damping while only slightly 
decreases the optomechanical resonant frequency by a 
few Hz (from the top dot to the two lower-right dots). 
Here we do not plot the optical resonances, which for the 
circulating power considered here remain stable. 

Before ending this section, we note that, when pendu- 
lum frequency is not neglected, there does exist a stable 
single optical spring regime where an increase in mechan- 
ical resonant frequency is associated with an increase in 
damping. But this requires that the optical frequency 
scales be lower than the pendulum frequency, which is 
not desirable in case of GW detectors. Such a regime 
was experimentally investigated by Schliesser et al. [2]| . 



III. INPUT-OUTPUT RELATION AND 
COMBINED NOISE SPECTRAL DENSITY 

In this section, we consider the sensitivity of stable 
DOS interferometers. Since the two carrier fields are dif- 
ferent in polarization and frequency, both fields can be 
measured independently and we effectively obtain two in- 
terferometers in one scheme, both sensing the same differ- 
ential mode [Eq. ([T])] of mirror motion. The input optical 
vacuum fluctuations associated with the two carriers, on 
the other hand, are independent. 

Throughout this paper, we will assume that GWs with 
amplitude h incident from right above the detector plane, 
with a polarization that maximizes the response of our 
L-shaped Michelson interferometer. In the following we 
will list the Heisenberg equations of motion in frequency 
domain [§, [l^, [ll|, |22| \23^ for the antisymmetric mode 
of motion of the arm cavity mirrors x and for the two 
measurement outputs y*-*' : 

+ (R^m + R%in)) X]+Lh + Cnoisc (13) 

=Yl^\n) sinC^^^ +Y^^\n) cosC*^^ 

+ [R^y^pm sinC^'^ + R'^^pin) cosC^i)]^, (14) 

y(2) ^Yl^\n) sinC'^'^^ +Y^^\n) cosC^^^ 

+ [R'y^pm sinC(2)-Fi?i|V(f^) cosC(^)]£. (15) 

Note that we have labeled all quantities with superscripts 
(1) and (2) for the first carrier and the second carrier, 
respectively. It was shown in Ref. [1^ that the interfer- 
ometer's output is only marginally influenced by seismic 
noise, thermal noise and radiation pressure noise intro- 
duced at the beam splitter, since the carrier light incident 
on the beam splitter is weak and the arm cavities prevent 
fluctuations from building up. The out-going sideband 
flelds at the dark port around the two different carrier 
flelds are split and each is sensed independently by a ho- 
modyne detection scheme, which measures a certain com- 
bination of amplitude and phase quadratures (described 
by the phases C^^^'^^^)- The operator F^') in Eq. 1^ de- 
scribes the radiation pressure forces which would act on 
flxed mirrors due to the incoming vacuum fields at the 
dark port. The operators Y^^^ in Eqs. (|14|) and (|15|) ac- 
count for the out-going fluctuations in the quadratures 
in case of fixed mirrors. (In language of Ref. [ll|, F^*' 

and Yj^*"* are free quantities). The operators ^noiso de- 
scribe the classical displacement noise of the differential 
mode. The quantity Rxx = —4:/{mfl'^) is the mechanical 
susceptibility of the differential mode, the susceptibili- 
ties R^'pp = —Kos' [cf. Eq. ^] correspond to the optical 

(i) 

spring constants, and Ry p are optical transfer functions 
from the differential mode to the out-going quadrature 
fields. 

According to Ref. [ll|, the free quantities F*-'^ and 
Y^*-*^ are related to incoming amplitude and phase quadra- 
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FIG. 4: Allowable regions for the optical resonances of the 
second carrier, given a circulating power of 8kW, 40 kW and 
80 kW, in order to stabilize a 800 kW first carrier that corre- 
sponds to the Advanced LIGO narrowband mode (cf. Tab. 
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For vacuum fluctuations, we have 

The optical transfer functions are given by [Tlj 
/fWeWm AW 
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(17) 
(18) 

(19) 
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2a (17 + ieW)2 - (AW)2 ' 
Finally, the classical noise operator ^noisc satisfies 

(Inoise(f^) ILsc(f^'))sym = TT <5(f7 - 17') L25;;'(f7) , (22) 

where S'f^ is the classical noise in the gravitational-wave 
strain. 

The outgoing sideband fields at the dark port around 
the two different carriers are detected independently via 
homodyne detection. We seek a linear combination of the 
two measurement output channels, y^^^ and y'^' given in 
Eqs. dUl) and (HH), 



(2) 



(23) 



which has optimal sensitivity to the GW strain h. 
Namely, one has to minimize the /i-referred noise spectral 
density by varying the two filter functions Ki(yi). The 
optimal solution can be found by a straightforward cal- 
culation described in detail for a general multi-channel 
interferometer in App. \^ 



IV. EXAMPLE CONFIGURATIONS 

In this section, different example configurations of a 
double optical spring interferometer are discussed. At 
first the considerations are restricted to quantum noise 
only in order to clarify two distinct regimes of our pro- 
posed scheme. Subsequently a configuration with real- 
istic classical noise budget is investigated and optimized 
for neutron-star-neutron-star binary inspirals. Finally, 
we lower the classical noise budget and explore whether 
DOS configurations can take full advantage of this im- 
provement. We are aiming at applying the DOS configu- 
ration as an upgrade candidate for the Advanced LIGO 
detector. 



A. Quantum noise examples 

Here we study the quantum noise spectrum of two spe- 
cial regimes of DOS: weak stabilization and annihilation. 

Weak stabilization. In this scenario, we use a rela- 
tively weak second carrier to stabilize a typical Advanced 
LIGO configuration. In Fig. [31 examples were already 
given, in which a second carrier of 8kW and 80 kW, 
respectively, has been used to stabilize the narrowband 
mode of Advanced LIGO (cf. Tab. |T|, while Fig. g] pro- 
vides the region of all possible second carriers, given cir- 
culating power of 8kW, 40 kW and 80 kW. 

In Fig. [5l we plot the noise spectrum of DOS interfer- 
ometers that correspond to the two stabilizing configura- 
tions in Fig. [31 namely with first carrier identical to the 
carrier of the Advanced LIGO narrowband mode, with 
resonant frequencies and powers of second carrier given 
in Fig. [31 and a phase readout quadrature associated with 
the second carrier (i.e., C'-^-' = 0). As we see from Fig. [SI 
the DOS noise spectra in both cases do not differ much 
from that of the Advanced LIGO narrowband mode. 

Annihilation. Now we turn to a different situation 
where both carriers have half the Advanced LIGO cir- 
culating power, namely 400 kW with the two detunings 
exactly opposite. In this case the two optical springs 
cancel each other and the total effective ponderomotive 
rigidity vanishes. For Fig.lHlthe two detection angles also 
have opposite signs while the absolute values of the two 
detunings (detection angles) agree with the detuning (de- 
tection angle) in case of the conventional Advanced LIGO 
configurations (cf. Tab.|T]). All other parameters are left 
unchanged. 
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detuning phase of SR cavity 




0^ n /I 
ZTT U.24/ 


PSR 


signal-recycling mirror reflectivity 




^093 


lo 


cavity half bandwidth 


27r f 5 Hz 


27r 15 Hz 




detection angle 


27r 0.347 


27r 0.45 


e 


effective half bandwidth 


27r 120 Hz 


27r 395 Hz 


A 


effective detuning 


27r 290 Hz 


2-K 411 Hz 



TABLE 1: Parameter values for single optical spring Advanced LIGO interferometer configurations used throughout the calcu- 
lations. The narrowband configuration is optimized for NS-NS binaries using the current Advanced LIGO noise budget. For the 
broadband operational mode, we allow 10% decrease in detectable distance for NS-NS binaries assuming the current Advanced 
LIGO noise budget, while maximizing the contribution to SNR from frequencies above 500 Hz 




f [Hz] 

FIG. 5: Noise spectrum of weak-stabilization DOS config- 
urations, with parameters corresponding to configuration in 
Fig. [3] Namely, the first carrier is identical to the carrier in the 
Advanced LIGO narrowband mode (cf. Tab. Two possible 
choices for the stabilizing second carrier are: (i) = 8 kW, 
e'^) = 27r 5Hz, A^^' = -27r 55Hz and (ii) P'^^ = 80 kW, 
e'^) = 2-K 60Hz, A<^' = -27r 60Hz. Phase quadrature detec- 
tion (CP'^ — 0) in both cases. 



The sensitivity does not change in the high frequency 
regime where shot noise is the limiting factor (cf. Fig.lH]). 
This is because, generally, the noise spectral density of 
the shot noise remains unchanged under reversing the 
sign of the detuning and the detection angle. With the 
optimal filter (cf. Fig. [7]), wc obtain the same shot noise 
as the Advanced LIGO configuration, as the total power 
is conserved. In the low frequency regime the sensitivity 
is slightly improved compared to the single optical spring 
interferometer. It is well-known that the quantum-noise 
limited sensitivity of single optical spring interferome- 
ters at frequencies below the optomechanical resonance 
is dramatically lower than the sensitivity of non-optical- 
spring interferometers because the strong restoring force 
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FIG. 6: Noise spectrum of DOS configurations with canceled 
optical springs, based on Advanced LIGO narrowband and 
broadband modes (cf. Tab. |l|. Detuning and detection angle 
of the first carrier are identical to that of Advanced LIGO 
narroband mode, while those of the second are opposite. The 
total power is equally divided into two parts [P^^^ = P*-^^ = 
400 kW). Other parameters are left unchanged. 

due to the single optical spring suppresses the response of 
the interferometer's differential mode to GWs. In a dou- 
ble optical spring interferometer the second carrier usu- 
ally results in a less rigid or completely canceled effective 
optical spring. This effect appears in the first example 
where only a weak second carrier is used (cf. Fig. [5]) and 
becomes more significant here in the case of canceled op- 
tical springs (cf. Fig. [6]). For intermediate frequencies the 
sensitivity becomes worse due to the absence of the op- 
tomechanical resonance gain. In this regime the output 
associated with the carrier having positive detuning usu- 
ally has more sensitivity than the other one. Therefore, 
as shown in Fig. [71 the filters make sure that only this 
output contributes to the total output. Note that when 
the two outputs are filtered appropriately, the total noise 
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FIG. 7: Example by means of the broadband configuration 
(cf. Fig.|6]): Squared absolute values of filter functions for the 
two outputs. These curves actually account for how much 
of the contributions is used - but they cannot show how the 
correlations among the two outputs have influence on the total 
output. The contributions correspond to measuring one of the 
outputs alone without filtering - but still in presence of both 



curve can in fact be below the single contributions. 

Since the optomechanical resonance disappears com- 
pletely, the noise spectral density in this configuration is 
equal to 



■^sql(^) 
4S',hot(f^) 



(24) 



Here S'shot(f^) denotes the shot-noise spectral density of 
a detuned single carrier interferometer with the same to- 
tal power as in the canceled optical spring configuration. 
The free-mass standard quantum limit (SQL) [25[ for de- 
tecting the gravitational-wave strain h with a Michclson 
interferometer with arm-cavitics is given by 



'SQL 



8h 



(25) 



B. 



Optimized configurations with Advanced LIGO 
classical noise buget 



In the following we optimize the double optical spring 
interferometer for neutron star - neutron star (NS-NS) 
binary inspirals. For such systems the last stable circular 
orbit gives an upper frequency limit of /max ~ 1570 Hz 
and seismic noise defines a lower bound of /min ~ 7 Hz, 
and the signal-to-noise ratio (SNR) in power achievable 
by optimal filtering, for a source at a given radius, is 
given by 



SNR^ cx 



/ 



-7/3 



Shif) 



d/. 



(26) 



Here Shif) denotes the single-sided noise spectral density 
of the interferometer. Note that this optimization strat- 
egy tend to focus more on the low frequency regime at 



the expense of the sensitivity at higher frequencies - due 
to the rather steep power law of f~'^^^. Since the radius 
of detectable range is proportional to the SNR at a fixed 
radius, and that the event rate is roughly the cube of the 
radius of detectable range, the event rate is proportional 
to SNR^ or 



Event Rate 



j-7/3 



d/ 



3/2 



(27) 



For the optimization we have taken the current Ad- 
vanced LIGO classical noise budget into account (as given 
in Bench H^): each contribution to the total classical 
noise budget, i.e. suspension thermal noise, seismic noise, 
thermal fluctuations in the coating and gravity gradient 
noise, arc presented in Figs.[5Hni 

It turns out that the second carrier not only helps to 
stabilize the interferometer but can also improve its sen- 
sitivity. Wc flrst assume both carriers to have the same 
SR mirror reflectivity, in this case the gain in NS-NS sen- 
sitivity is maximized when the two optical springs totally 
cancel each other (Fig. [5]) at an equal power distribution 
in the two carrier fields. Fig. [8] also shows the contribu- 
tion of each carrier to the total noise spectral density. 

If each carrier senses a different refiectivity of the SR 
mirror the interferometer's sensitivity can be improved 
further. With this additional degree of freedom the sit- 
uation of canceled optical springs is no longer the op- 
timal choice for NS-NS binaries. Each carrier can be 
optimized for a different frequency regime such that they 
complement each other. This is illustrated for an equal 
distribution of the total power (400 kW for each carrier) 
in Fig. [5] where one carrier ensures good sensitivity in 
the low frequency regime while the other one gives the 
main contribution at frequencies above the optomechan- 
ical resonance as shown by the contribution curves. The 
noise spectral density is close to the classical noise level 
in low frequencies, and an improvement in event rate of 
35.5% can be achieved for NS-NS binaries for which the 
optimization was performed. As a more general opti- 
mization is performed, allowing two carriers to have dif- 
ferent powers summing up to 800 kW, it turns out that 
= 500 kW, P(2) = 300 kW achieves a slightly better 
event rate improvement of 36.3%. Results for different 
power distributions are given in Tab. [Ill 

The sensitivity of all optimized configurations is basi- 
cally improved at low frequencies at the expense of the 
high frequencies. Therefore it might be necessary to carry 
out an additional broadband optimization in order to 
achieve a better sensitivity in the high frequency regime. 
This can be accomplished by first picking out all config- 
urations obeying an event rate which is at least a certain 
fraction of the optimal event rate for NS-NS binary sys- 
tems (cf. Eq. (HH)). In a second step these configurations 
are explored in the high frequency regime by considering 
a smaller frequency integration interval [500 Hz, 1570 Hz] 
and selecting the optimal signal-to-noisc ratio on this in- 
terval. For the (F^^^ = P'^' = 400 kW)-configuration as 
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P(i) in kW P(^^ in kW A'^' in Hz A'^' in Hz e^^^ in Hz e'^^ in Hz C'^^ in radian C'^' in radian Improvement 





800 





2-K 290 




2n 120 




27r 0.347 






'3 

e 


750 


50 


2n 190 


-27r 70 


27r 50 


27r 20 


27r 0.125 


27r 0.075 


7.3% 


ical 


700 


100 


27r 190 


-27r 50 


27r 40 


27r 30 


27r 0.0625 


27r 0.475 


22.6% 


ra 


600 


200 


2n 180 


-2n 40 


2-K 30 


2k- 50 


2-K 0.0125 


27r 0.3875 


34.4% 


"3 


500 


300 


27r 150 


-27r 50 


27r 45 


27r 60 


27r 0.0125 


27r 0.375 


36.3% 




450 


350 


2-rr 160 


-2n 10 


2-K 30 


27r 40 


27r 


27r 0.3 


35.5% 




400 


400 


2n 140 


-27r 20 


27r 20 


2-K 55 


2-K 0.01 


2-K 0.31 


35.5% 




800 





2ir 170 


- 


27r 10 


- 


27r 0.263 


- 


- 


'3 

a 


700 


100 


2tt 165 


-2tt 75 


2tt 5 


27r 25 


2-K 0.48 


27r 0.02 


32.1% 


"3 


600 


200 


27r 150 


-27r 20 


27r 5 


2tv 45 


27r 0.01 


2% 0.36 


52.2% 




500 


300 


2tt 140 


-27r 20 


2tt 5 


27r 45 


27r 0.03 


27r 0.33 


74.3% 




450 


350 


271 135 


-2tt 25 


27r 5 


2tv 50 


2tv 0.01 


27r 0.33 


83.4% 




400 


400 


2tt 127 


-27r 11 


27r 4 


27r 46 


27r 0.01 


27r 0.3 


110% 



TABLE II: Parameters for double optical spring scheme optimized for NS-NS binary systems. The total power is fixed to 
800 kW. The last column gives the improvement in event rate for our proposed scheme compared to an optimized Advanced 
LIGO configuration provided in the first row. For the upper part we adopted the current Advanced LIGO noise budget and 
for the lower part the gravity gradient noise and the suspension thermal noise are reduced by a factor of ten and the coating 
thermal by a factor of three. 




f [Hz] 

FIG. 8: Comparison of single optical spring and double carrier 
with canceled optical spring Advanced LIGO configuration: 
the same refiectivity of the SR mirror is assumed for both 
carriers. We obtain 15% improvement in event rate for DOS 
interferometer using the following parameters: P^^^ = P'-'^' — 
400 kW, p'il^ = p^ily^ = 0.87, = -0(2' = 27r 0.233, 
C'^' = 2-K 0.433 and C*^' = 0. 




f [Hz] 

FIG. 9: Comparison of single and double optical spring Ad- 
vanced LIGO configuration: different SR mirror reflectivities 
for first and second carrier are allowed. Parameters are given 
in detail in last row of first block in Tab.lIIlfor the narrowband 
configuration. For the broadband we have: A'^' = 2-k 130 Hz, 
e'l' = 27r 5 Hz, C^^' = Stt 0.01, A'^' = -2k: 230 Hz, 
e'2' = 2-K 155 Hz and C*^^ = 27r 0.02. We obtain 35.5% 
improvement in event rate in case of the narrowband config- 
uration and 16.4% for the broadband configuration. 



an example, we can achieve a sensitivity comparable to 
Advanced LIGO on the [500 Hz, 1570 Hz] frequency band 
while maintaining an improvement in the event rate (inte- 
grating again from /min) of 16.4% compared to Advanced 
LIGO (cf. FiglH). 

After adding up all classical noise contribution shown 
in Figs.[8]and[9]it turns out that the noise spectral den- 
sity of the double optical spring configuration - contrary 
to the single optical spring in Advanced LIGO - almost 
follows the borderline set by the classical noise in the low 



frequency regime. Let us write 

|.150Hz , |.150Hz y-7/3 



Qcl 



Ccl 



+ 5? ' 



(28) 



where S'f^, S'^ and Sh are weighted averages of classical, 
quantum, and total noise spectrum in the low-frequency 
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band, respectively. For the P^^^ = P^^^ = 400 kW con- 
figuration (cf. Fig.[ni) we obtain 

?7«0.81. (29) 

This indicates that, at low frequencies (i.e., below 
150 Hz), the quantum noise is already small fraction of 
the total noise; improving quantum noise further does not 
significantly improve sensitivity. Qualitatively, starting 
with a level of r/ ~ 0.81, further lowering S'^ by a factor 
of 2 only improves r] from 0.81 to 0.88, which yields an 
14% increase in event rate. 

C. Optimized configurations with classical noise 
budget beyond Advanced LIGO 

It is likely that technical improvements will reduce the 
classical noise floor in the future. In order to explore 
the potential of our proposed configuration to increase 
the quantum limited sensitivity by reshaping the noise 
curves in an optimal way, we analyze the performance 
for a reduced classical noise budget. For instance, the 
gravity gradient noise is a limiting factor at lower fre- 
quencies. As suggested in Ref. [l^j this effect can be re- 
moved from the recorded data by performing an indepen- 
dent measurement of the ground's density fluctuations 
near each test- mass. We assume it to be 1/10 (in am- 
plitude) the current estimation for Advanced LIGO [1^. 
Another limiting factor is given by the thermal noise in 
the suspension system and mirrors. We assume that the 
suspension thermal noise can be lowered by a factor of 
10 in amplitude, while internal thermal noise of mirrors 
can be lowered by a factor of 3 in amplitude [l^ . Such 
improvements may possibly be realized by (i) optimiz- 
ing the design of mirror coating structure and suspen- 
sion wires, (ii) improving mechanical quality factors of 
mirror coating, substrate and suspension materials, and 
(iii) applying cryogenic techniques [13, [H, HO, [3l|, HI] . 

Now we optimize the usual single optical spring con- 
figuration as well as our proposed double optical spring 
layout for this modified noise budget. Only now the po- 
tential of the two carrier interferometer becomes appar- 
ent: For the new classical noise budget, the single optical 
spring configuration can improve the event rate for NS- 
NS binary inspirals only by 61%. But the DOS configura- 
tion can do 238%. This corresponds to an improvement 
in event rate of the DOS compared to the optimized (with 
respect to the new classical noise budget) single optical 
spring configuration of 110% (cf. Tab. HI)) . 

Such a situation is presented in Fig. 1101 The big gap 
between the single optical spring total noise and the clas- 
sical noise budget verifies that the single optical spring 
Advanced LIGO configuration has not been limited by 
classical noise at low frequencies. While this gap can 
be partially filled by the double optical spring configura- 
tion, there is still room for further improvement. For the 
configuration used in Fig. [TO] we evaluate 

77 w 0.3. (30) 




f [Hz] 

FIG. 10: Compared to Fig.[H] suspension thermal noise and 
gravity gradient noise are lowered by a factor of 10 and coat- 
ing thermal noise by a factor of 3. DOS configuration as well 
as single optical spring configuration both optimized with re- 
spect to the new noise budget. An improvement in event 
rate of 238% can be achieved by DOS configurations, com- 
pared with 61% of single optical spring configurations. The 
following parameters were used: P'^' = P'^-* = 400 kW, 
A'l' = 2tv 130 Hz, A(2) = -2tv 10 Hz, e^^^ = 2tv 5 Hz, 
e(2) = 2-K 45 Hz, C'^- = 27r 0.46, C'^^ = 2% 0.3. For single 
optical spring optimization, parameters are given in the first 
row of the second block in Tab. [Hi 

In this case, lowering the quantum noise by a factor of 2 
improves the event rate by 69%. One possibility of fur- 
ther improvement would be to inject even more than two 
carriers, and combine the corresponding output channels 
(cf. App.E|. 

V. CONCLUSION 

While the concept of a stable Double Optical Spring 
(DOS) has motivated experiments such as the one al- 
ready carried out by Corbitt et al. [l^], in this paper, 
we have theoretically investigated the benefit of a Dou- 
ble Optical Spring (DOS) configuration for second gen- 
eration gravitational-wave detectors, in particular in a 
follow-up experiment to Advanced LIGO, possibly in 
combination with other existing schemes, e.g., the local 
readout scheme proposed in Re f. II6II1 as well as the in- 
jection of squeezed vacuum [H, |34|. l35l. [36l. ISTf. 

In the DOS configuration, a second laser beam is in- 
jected into signal-recycling interferometers at the bright 
port which is, as the first carrier, resonant in the arm- 
cavities and is also detuned in the signal-recycling cav- 
ity. The two outputs are optimally filtered and combined. 
By choosing appropriate detunings of the signal-recycling 
cavity and homodyne detection angles, it is possible to 
achieve a stable double optical spring while additionally 
improving the sensitivity. 

Taking into account the current classical noise bud- 
get estimation for the Advanced LIGO detector, as well 
as constraints on optical power, we have performed an 
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optimization of our double-carrier scheme towards the 
detection of compact binary inspirals specifying to neu- 
tron stars. The DOS allows a 36% improvement in 
event rate, and we have shown that further improvement 
in event rate will largely be limited by classical noise. 
When considering a more optimistic classical noise bud- 
get, DOS interferometers are much more capable than 
singlc-optical-spring interferometers in taking advantage 
of this improvement: compared with 61% improvement 
in event rate achievable by single-optical-spring config- 
urations, the DOS allows 238%. Nevertheless the re- 
duced classical noise level leaves further room beyond 
DOS, which can be exploited by employing more than 
two carrier fields. 



where i7 is a vector with 2n + l entries which account for 
2n quadrature operators and one operator modeling the 
classical noise. Here Hi describes the noise transfer func- 
tions from the noise channels v into the output channels 
and Si accounts for the signal-transfer functions. The 
combined output is given by 



y(r!) = ^if,(f^)y«(0), (A3) 

1=1 

and one has to identify n optimal filter functions which 
minimize the signal-referred noise spectral density of y: 
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APPENDIX A: OPTIMAL OUTPUT CHANNEL 
FOR INTERFEROMETERS WITH MULTIPLE 
CARRIERS 

Here we analyze the optimal way of combining the out- 
put channels obtained from homodyne detections at the 
dark port of a multi-carrier interferometer. The general 
equation of motion for the mirror position in case of n 
carriers reads [cf. Eq. ([T^ ]: 



Shin) 



with 



and 



2n+l 
s,k=l 



(A4) 



(A5) 



(A6) 



Here S^^,y^ is the cross spectral density between i/g and 
i/fc, with 



(A7) 



n 



(Al) 

The output corresponding to each of the n carrier fields 
is given by [cf. Eqs. (fH)) and ([T5])]: 

y^'\n)=Yl'\n) sinC^'^ +Y2^'\n) cosC^'^ 

+ sinC^'^ + Ri;)p{n) cosC(')]x 



=fij u + Sih 



(A2) 



The inverse of the largest eigenvalue of the n x n- matrix 



M = • S. 



(A8) 



provides the resulting minimum noise and the corre- 
sponding eigenvector gives the n optimal filter function 



[1] D. Shoemaker, and The LIGO Scientific Collaboration, [2] L. D. Fiore, and VIRGO collaboration, 
Nucl. Inst, and Meth. A 517, 154 (2004). 



11 



Class. Quant. Grav. 19, 1421 (2002). 
[3] B. Willke, et al, Class. Quantum Grav. 19, 1377 (2002). 
[4] M. Ando, and TAMA collaboration, Phys. Rev. Lett. 86, 

3950 (2001). 

[5] Advanced ligo reference design, report ligo-060056-07-m 
(2006), http://www.ligo.caltecli.edu/docs/M/M060056- 
07/M060056-07.pdf. 

[6] B. Meers, Phys. Rev. D 38, 2317 (1988). 

[7] J. Vinet, B. Meers, C. Man, and A. Brillet, Phys. Rev. D 
38, 433 (1988). 

[8] G. Heinzel, K. A. Strain, J. Mizuno, K. D. Skeldon, 
B. Willke, W. Winkler, R. Schilling, A. Riidiger, and 
K. Danzmann, Phys. Rev. Lett. 81, 5493 (1998). 

[9] A. Buonanno and Y. Chen, Phys. Rev. D 64, 042006 

(2001) . 

[10] A. Buonanno and Y. Chen, Phys. Rev. D 65, 042001 

(2002) . 

[11] A. Buonanno and Y. Chen, Phys. Rev. D 67, 062002 

(2003) . 

[12] K. Somiya, P. Beyersdorf, K. Aral, S. Sato, S. Kawamura, 

O. Miyakawa, F. Kawazoe, S. Sakata, A. Sekido, and 

N. Mio, Appl. Opt. 44, 3179 (2005). 
[13] O. Miyakawa, et al, Phys. Rev. D 74, 022001 (2006). 
[14] V. B. Braginsky, Y. I. Vorontsov, and F. Y. Khahh, Sov. 

Phys. JETP 46, 705 (1977). 
[15] V. B. Braginsky, M. L. Gorodetsky, and F. Y. Khahli, 

Phys. Lett. A 232, 340 (1997). 
[16] H. Rehbein, H. Mueller-Ebhardt, K. Somiya, C. Li, 

R. Schnabel, K. Danzmann, and Y. Chen, Phys. Rev. D 

76, 062002 (2007). 
[17] E. Gustafson, D. Shoemaker, K. Strain, 

and R. Weiss, LIGO-T990080-00-D (1999), 

http://www.ligo.caltech.edU/docs/T/T990080-00.pdf. 
[18] T. Corbitt, Y. Chen, F. Khahh, D. Ottaway, S. Vy- 

atchanin, S. Whitcomb, and N. Mavalvala, Phys. Rev. A 

73, 023801 (2006). 
[19] H. Mueller-Ebhardt, H. Rehbein, T. Corbitt, C. Wipf, 

N. Mavalvala, R. Schnabel, K. Danzmann, and Y. Chen, 

in preparation (2008). 
[20] T. Corbitt, Y. Chen, E. Innerhofer, H. Mueller-Ebhardt, 

D. Ottaway, H. Rehbein, D. Sigg, S. Whitcomb, C. Wipf, 

and N. Mavalvala, Phys. Rev. Lett. 98, 150802 (2007). 



[21] A. Schliesser, P. Del'Haye, N. Nooshi, K. J. Vahala, and 
T. J. Kippenberg, Phys. Rev. Lett. 97, 243905 (2006). 

[22] H. J. Kimble, Y. Levin, A. B. Matsko, K. S. Thorne, and 
S. P. Vyatchanin, Phys. Rev. D 65, 022002 (2001). 

[23] T. Corbitt, Y. Chen, and N. Mavalvala, Phys. Rev. A 
72, 013818 (2005). 

[24] J. Harms, R. Schnabel, and K. Danzmann, Phys. Rev. D 
70, 102001 (2004). 

[25] V. B. Braginsky and F. Y. Khalili, Quantum Measure- 
ment (Cambridge University Press, 1999). 

[26] http://www.ligo.mit.edu/bench/bench.html. 

[27] S. A. Hughes and K. S. Thorne, Phys. Rev. D 58, 122002 
(1998). 

[28] R. Adhikari, private communication (2007). 

[29] A. Bunkowski, O. Burmeister, D. Friedrich, K. Danz- 
mann, and R. Schnabel, Classical and Quantum Gravity 
23, 7297 (2006). 

[30] M. Bondarescu and K. S. Thorne, Phys. Rev. D 74, 
082003 (2006). 

[31] T. Uchiyama, T. Tomaru, M. E. Tobar, D. Tatsumi, 
S. Miyoki, M. Ohashi, K. Kuroda, T. Suzuki, N. Sato, 
T. Haruyama, et al., Phys. Lett. A 261, 5 (1999). 

[32] T. Uchiyama, D. Tatsumi, T. Tomaru, M. E. To- 
bar, K. Kuroda, T. Suzuki, N. Sato, A. Yamamoto, 
T. Haruyama, and T. Shintomi, Phys. Lett. A 242, 211 
(1998). ' 

[33] K. McKenzie, N. Crosse, W. P. Bowen, S. E. Whit- 
comb, M. B. Gray, D. E. McClelland, and P. K. Lam, 
Phys. Rev. Lett. 93, 161105 (2004). 

[34] H. Vahlbruch, S. Chelkowski, B. Hage, A. Franzen, 
K. Danzmann, and R. Schnabel, Phys. Rev. Lett. 97, 
011101 (2006). 

[35] J. Harms, Y. Chen, S. Chelkowski, A. Franzen, 
H. Vahlbruch, K. Danzmann, and R. Schnabel, 
Phys. Rev. D 68, 042001 (2003). 

[36] A. Buonanno and Y. Chen, Phys. Rev. D 69, 102004 
(2004). 

[37] H. Vahlbruch, S. Chelkowski, B. Hage, A. Franzen, 
K. Danzmann, and R. Schnabel, Phys. Rev. Lett. 95, 
211102 (2005). 



